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Summary.  In  this  report  we  consider  the  number  N  of  uperossings  of  a  level  u  by  a 
stationary  normal  process  |(t)  in  0  ^  t  ^  T.  A  formula  is  obtained  for  the  factorial 
moment  ■  £{N(N- 1 ) . . . (N-k+1 )  )  of  any  desired  order  k.  The  main  condition  assumed 
in  the  derivation  is  that  P(t)  have,  with  probability  one,  a  continuous  sample  deri¬ 
vative  ?’(f)  in  the  interva 1  [0 ,T] .  This  condition  involves  no  real  restriction  since 
an  example  shows  that  even  a  slight  relaxation  of  it  causes  all  moments  of  order 
greater  than  one  to  become  infin'te.  The  moments  of  the  number  of  down  crossings  or 
total  number  of  crossings  can  be  obtained  analogously. 

1.  Introduction.  The  problem  of  obtaining  the  mean  number  of  crossings,  (or  equi¬ 
valently  uperossings)  of  a  given  level,  by  a  stationary  normal  process  in  a  given  time 
has  received  a  good  deal  of  attention  in  the  literature.  In  fact,  a  complete  solution 
to  this  problem  has  now  been  given  by  Ylvisaker  [8].  However,  moments  of  order  greate 
than  me  of  the  number  of  crossings  of  a  level  have  received  less  attention.  The 
variance  was  obtained  by  Steinberg  et  al  [6],  using  somewhat  heuristic  arguments. 
Rozanov  and  Volkonski  [7]  point  out  in  a  footnote  that  the  formula  given  in  [6]  for 
the  variance  is  valid  under  certain  precise  conditions,  of  which  the  main  one  is 
that  the  covariance  function  of  the  process  have  a  finite  sixth  derivative  at  the 
origin.  Finally  in  this  connection,  the  variance  has  been  obtained  by  Leadbetter  and 
Cryer  [4]  under  conditions  which  assume  just  a  little  more  than  the  existence  of  a 
second  derivative  of  the  covariance  function. 

Ihere  is  virtually  no  literature  available  in  connection  with  moments  of  the 
number  of  crossings  of  a  level,  of  higher  order  than  the  second.  (  A  partial 
result  is  indicated  by  Ivanov  at  the  end  of  his  paper  [3]  ).  It  will  be  our 


purpose  here  to  obtain  explicit  expressions  for  such  moments.  This  will  be  done  for 
uperossings,  in  terms  of  factorial  moments  of  arbitrary  orders  and  under  conditions 
which  are  very  close  to  the  necessary  ones.  Corresponding  formulae  for  moments  of  all 
orders  for  the  down  crossings,  or  total  number  of  crossings,  follow  similarly. 

2.  Moments  of  the  number  of  uperossings.  We  shall,  throughout,  consider  a  real 
valued  stationary  normal  process  (  5(t):  0/t<T)  having  (for  convenience)  zero 
mean,  spectrum  F(X)  possessing  an  absolutely  continuous  component,  and  covariance 
function  r(r)  =  /  e1'^‘dF(X).  We  shall  further  assume  that  |(t)  has,  with  probabi  Uty 

one,  a  continuous  sample  derivative  |*(t)  on  the  interval  [0,T].  Sufficient  conditions 
for  this  latter  property  in  terms  of  the  behaviour  of  the  covariance  function,  are 
well  known.  Write  N  for  the  number  of  uperossings  of  the  level  u  by  |(t)  in 
0  <  t  <  T;  that  is  N  is  the  number  of  points  t  in  that  interval  for  which 
P,(t)  =  u,  P'(t)  >  0.  Then  the  following  result  holds. 

Theorem. 

If  (f,(t):  0  <  t  <  T]  is  a  normal  stationary  process,  as  described,  possessing, 

with  probability  one,  a  continuous  sample  derivative,  and  k  is  any  positive  integer, 
then 

Mk  =  £{N(N-1). . . (N-k+1) ) 

T  *• 

(1)  =  /.../  dtj . . *dtR  /.../  y1...yR  Pt(u,£)  dy^ . . . dyk 

o  o  - 

in  which  Pt(u,i>  ■=  pfc (u  ,  .  .  .u.yj  ,  .  .  .yk>  ,  p(.(x1 , .  .  .xk?yl .  .  .yR)  denoting  the  joint 

density  for  the  random  variables  £ (t^) . . . 5 (tk)  ,  | 1  (t^ ) . . . £ ' (tk) .  We  note  here  that 
it  follows  from  the  appendix  that,  when  all  t^  are  different,  this  is  the  density 
corresponding  to  a  non  singular  joint  distribution  since  F(X)  is  assumed  to  have  an 
absolutely  continuous  component. 
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Before  proceeding  to  the  proof  we  note  that  the  theorem  can  easily  be  modified 
to  refer  to  "downcrossings"  or  the  total,  number  of  crossings  of  the  level  u  in  time 
T.  The  discussion  will  be  given  here  in  terms  of  uperossings,  however. 

The  following  proof  is  divided  into  two  parts  A  and  B.  In  Part  A  it  is  shown 
that  does  not  exceed  the  expression  on  the  right  hand  side  of  (l),  whereas  in 
Part  B  the  reverse  inequality  is  proved.  The  techniques  are  straightforward  but 
quite  different  in  each  part.  It  is  a  perhaps  somewhat  surprising  feature,  however, 
that  in  both  parts  use  can  be  made  of  Fatou's  Lemma  to  give  the  essential  inversions 
of  limiting  operations  with  integrations,  in  order  that  inequalities  in  the  desired 
(opposite)  directions  may  be  obtained. 


Proof  of  the  Theorem,  Part  A. 

Write  5(t)  =  ^(t,u))  to  exhibit  explicit  dependence  on  the  "sample  point"  (a  e  £2. 
Let  S  denote  the  set  of  all  oo  such  that  the  equation  ?(t)  =  u  has  at  most  a  finite 
number  of  roots  t  in  the  interval  I  ■  [0,T],  while  further  |(  0;  4  u  +  £(T)  and 
*'(t)  4  0  whenever  f(t)  =  u.  According  to  Bulinskaya  f 2] ,  Theorem  1  we  then  have 


(2) 


Write  now  N 


(k) 


P(S)  =  1 

N(N- 1 ) . . . (N-k+1 )  for  k  =  1,2...  and  define  the  functions 


b  (x),  a(x)  by 
n 


and 


•  n(x)  =  n  |x  I  <  l/(2n) 


0  otherwise 


o(x)  =  x 


x  >  0 


=  0 


otherwise 


Let  D(e)  denote  the  domain  in  the  k-dimensional  space  R  with  coordinates 


t^.-t^  defined  by  the  inequalities 
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0  t .  <  T  for  i  =  1 . . .k 

i 

1 1  -t  .  |  >  r  for  i  4  j . 
i  J 

Define  also  the  random  variable  J^(n,c,oj)  by  the  relation 

k 

(3)  J.  (n,e,u>)  =  /.  . ./  7  ~  [ £  ( t  ) -u]  o  [ £ '  ( t  )  )  dt  .  .  .  dt 

k  D(r)  i=l  “  1  l  i  K 

We  shall  now  proceed  to  prove  that 

(4)  -'{N(k))_  lim  lim  t'(Jk(n,e,oo))  . 

e  — >  o  n  — *>  *• 

In  order  to  prove  the  validity  of  (4)  we  define  a  subset  S(h)  of  S  consisting 
of  all  ua  c  S  for  which  the  following  two  conditions  are  satisfied 

(a)  The  distance  between  any  *-wo  upcrossings  of  :(t)  with  the  level  u  in  I 
is  g  eater  than  2h, 

(b)  For  any  zero  t  =  t  of  the  derivative  :’(t)  in  I,  we  have  |:(t  )  -  u!  >  h . 

o  o 

According  to  the  definition  of  S  every  co  c  S  must  also  belong  to  S(h)  for  some 
h  >  0.  For  it  is  obvious  that  property  (a)  is  satisfied  for  :(t,co)  if  h  is 
suf f i ci ent ly  small  and  if  (b)  were  not  satisfied  we  could  find  a  sequence  of  points 
t^  r  I  for  which 

£,(t.)  =  0  |*(t  )  -  u |  <  1/i  . 

But  such  a  sequence  (t-1  must  have  a  limit  point  t  e  I  and  the  continuity  of  -  and 

e'  show  that  :'(t  )  =  0,  £(t  )  =  u  contradicting  the  fact  that  j)  e  S.  Hence  we 
o  o 

thus  see  that 

(5)  S(h)  f  S  as  1  |  0 
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Take  now  any  fixed  x  c  S(h),  and  let  t  =  r^...r^  be  all  the  upcrossings  of 

the  corresponding  :(t)  =  :(t,co)  in  I.  Consider  the  k-dimensiona  1  interval  I  in 

,  k 

the  space  R  ,  and  let  A.  .denote  the  point  in  I  with  coordinates 

Jl’‘"Jk 


k 

where  each  j  may  assume  the  values  1,2...N.  Clearly  there  are  N  different  points 

(k) 

A,  and  among  these  there  are  exactly  N  points  A'  such  that  no  two  of  the  j  ^  are 
equal.  Since  o>  e  S(h),  these  points  A'  will  all  be  situated  in  the  domain  D(2h), 
while  the  remaining  Nk  -  points  A  will  fall  outside  D(2h),  and  even  outside 

D(r) ,  for  any  e  >  0. 

Considering,  still  the  same  fixed  co  c  S(!t)  we  now  take  n  and  c  such  that 

„  -1 

0  <  n  <  e  <  h 

and  consider  the  integral  J^(n,e,co)  defined  by  (3).  The  contribution  to  J^(n,r,i.) 

arising  from  small  disjoint  k  dimensional  blocks  about  each  point  A'  is  easilv 

seen  to  be  just  for  all  sufficiently  large  n  (i.o.  a  unit  contribution  from 

each  such  block).  The  contribution  from  the  remaining  region  is  zero  for  all 

sufficiently  large  n.  (This  can  be  seen  clearly  from  a  picture  by  taking  k=2 

and  writing  down  the  integrals  involved).  Hence  for  any  fixed  c  <  h,  we  can  always 

find  n  so  large  that,  for  all  n  >  n  we  have 
o  0 

Jk(n,c,u>)  = 

and  hence  also 

N(k)  =  Urn  J R ( n , e  , to) 

n  — 

Since  this  holds  for  any  c  •  h,  while  the  first  member  is  independent  of  r,  it  follows 
that 
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(6) 


lim 


J  ^  ( n ,  £  ,  co) 


N 


00 


lim 


e  — *>  u  n  — >  «• 

for  any  fixed  os  e  S (h ) .  But  h  can  be  chosed  auitrarily  small  and  since 
S(h)  'f  S  as  h  ^  0  it  follows  that  (6)  holds  for  any  cd  e  S,  i.e.  with  probability 
one.  Finally  an  application  of  Fatou's  Lemma  to  the  e  and  n-limits  yields  the- 
result  (4).  Thus  from  (4)  we  obtain 


(7)  t’{N(k)]< 


lim  lim  /. . . /dt 
e  — >o  n  — D(e) 


dtR[nk  /.  .  ./dxr 
uJL 

pn 


dxk:.  yk  pt(x1...xk,y1...yk)dy1...dyk] 


The  entire  expression  in  square  brackets  on  the  right  hand  side  of  (7)  clearly 

m 

converges  to  /.../ y^.-.y^  Pt(u,^)  dy^...dy^.  Further,  it  can  be  readily  shown  that 
o  — 

this  expression  is  bounded  fox'  ail  t^..t^  in  the  region  D(e)  (using  the  fact  that 
the  determinant  of  the  covariance  matrix  of  |(t^) .  . . |(t^)  ?'(t^)..  is  bounded 

away  from  zero).  Hence  by  dominated  convergence 


(8) 


£(N(k))  <  lim 

c — >o 


/  /  dt  .dt  J.../y  y  p  <u,Z)dy  ...dy 

D(e )  o 


Finally  by  monotone  convergence  it  follows  that 


(9) 


e(N(k) )  < 


■/dt, 


.dt. 


/••■/yl-..yk  Pt(u,l)dy1. . .dyk 
o  —  " 


Part  B. 

In  order  to  prove  the  reverse  inequality  to  (9)  we  adopt  a  different  procedure 
(due  to  Ylvisaker  [8))  for  counting  the  number  of  upcrossings  by  *(t)  in 
0  <  t  <  T.  First,  however,  we  note  that  if  X  ,  i  ■  1,2...  are  each  either  zero 
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m 


or  one,  and  M  ■  E  X  ,  then,  for  any  integer  k  <  m, 
i-1 


(10) 


M 


(k) 


M(M-l). . .(M-kfl)  -  E'X  ...X 

Xk 


where  E*  denotes  summation  over  all  possible  ordered  sets  of  distinct  integers 

i  ...i  .  For  M  is  just  the  number  of  non  zero  Xi  ,  and  the  right  hand  side  of  (10) 

1  k 

therefore  represents  the  number  of  ordered  sets  of  distinct  integers  i^.  .  .ik 
such  that  each  corresponding  XA  is  non  zero,  taken  out  of  a  total  of  M  possible 
integers  i  for  which  Xt  +  0.  Bit  this  number  is  simply  M(M-l) . . . (M-k+1)  as  required 
Write  now  ^  =  |(Ti/fn),  i  *  0,  1...2  ,  n  =  1,2^..  .  Let  X^  *  1  if 
| i  <  u  <  p 1+i»  and  Xi  *  0  otherwise.  Then  if  ■  E  ^  we  have  Nr  f  N  a.s.  . 

(A  detailed  proof  of  this  latter  statement  is  given  by  Ylvisaker  [8]).  Hence  by 
monotone  convergence, 


(ID 


(k) 

«(Nn  ) 


00 

^{N  }  as  n  - > 


Now  from  (10)  we  have  with  m  *  2  , 


(12) 


(k) 

g(Nn  ) 


E'  p(V.  -  X 

ir..ik  1  2 


X  -  1} 

k 


We  note  that  no  terms  for  which  |  i  - i  |  =  l  for  any  r,s  appear  since  we  cannot  have 

IT  S 

Xi  -  Xt  +1  -  1.  Write  tj1  -  2n({i+1-'i)/T.  Then 
r  r 


»(X  -  X .  X  *  1}  -  P(u-2  V)  <  <  u,  r  -  1  •  •  •  k ) 

ll  i2  *k  r  r 


/.../  dyv..dyv 


r  Vij 


(x1...xk,y1...yk)dx1...dxk, 


where 


u-2""^  u-2_nTyk 
p  Is  the  joint  density  for  the  distribution  of  ?  .  .  ,  T,  ■  ••rh 

n’ij  ll  k  ll  I 
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(That  this  distribution  is  non  singular  follows  from  the  Appendix.)  By  a  change 
of  the  x-variables  in  this  expression  we  thus  obtain  from  (12) 


(13)  <?fN<k)J-  2'knTkZ'/..Jdyi...dyk  f° .  .  f°  p  (u+2"1^, . .  •u+2‘nT3r  ^  ,  .  .  .yk) 

_yl  -yk  j 


dx, .  .  .  dx. 
1  k 


Write  now  ?n  £(*l .  .  .x^ .  .  .yfc)  -  Pn  ,  <*!  •  •  V>V  •  -Vw) 


for  all  t.-.t.  such  that  t  lies  in  the  interval  (1  T/2n,  (1  +l)T/2n)  for  each  r. 
Ik  r  r  r 

Then  (13)  may  be  rewritten  as 


o  o 


(14)  /.  .  ./ dtj .  .  .dtk  /.  •  ./dyr  .  .dyk  /.  .  ./  ?nt  (u+2 '"TXj^  , .  .  .u+2  "l^.yj.  .  .yk)dxr  .  .dx^ 


■yl  -yk 


where  R  is  the  subset  of  I'  for  which  no  two  of  t,...t,  are  contained  in  the  same 
o  lk 

or  adjacent  intervals  of  the  form  (r  T/2n,(r+l)  T/2n).  (See  the  remark  following 
Equation  (12)). 

Let  now  (t,...t,)  be  a  fixed  point  in  R  .  Then  ¥  is  a  2k-dimensional  normal 
1  k  o  nt 

density  function.  Suppose  that  i /2 n  <  t  <  ( i 1 ) T / 2 n ,  r»l...k.  Then  corresponding 

to  the  point  t^  we  have  the  random  variables  Ki^T/2  ),  T](i^f/2‘),  yielding  the 

following  typical  members  of  the  covariance  matrix  for  T  for  example 

nt 


var  (| .  )  -  r(0), 


writing  *  for  |(iT/2 ) 


c°v  •  r 


l  2 


writing  for  r(iT/2  ),  p  -  i ^ - i 2 


cov  (P  ,  v  )  ■  2n(l-r.)/T 
11 


c°v  (p^,  niz)  -  2n(Vl  -  rp)/T 

var  tj  -  2^n+1(l-r  )/T2 
1 

cov  (-l^.lj  )  »  -22n[r^i  '  2rp  + 
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For  the  fixed  t^,t2  considered  i^.i^.p  depend  on  n.  It  is  an  easy  exercise  to  show 

that  if  t  =  1 2  ■  t  ^  the  above  elements  converge  (in  the  order  given),  as  n  — ■>  •»,  to 

r(0),  r(i),  0,  r'(T),  -r"(0),  -r"(x),  respectively.  Similar  conclusions  hold 
for  the  elements  corresponding  to  any  pair  t^,  t  ^  .  But  this  means  that  the  integrand 

in  (14)  must  converge  to  p  (u,^)  as  n  — ■>  «•  and  hence,  by  Fatou  s  Lemma 

(k)  T 

(15)  ’  }  >  /.  ■  ./dt^  .  .dtk  /.  .  ./  Pt(u,^)dy1.  .  .dyR 

o  - 

Combining  (9)  and  (15)  we  obtain  the  desired  equality  and  hence  the  truth  of  the 
theorem  follows. 

3.  A  case  when  *  +  m.  Formula  1  was  obtained  under  the  condition  that  F(t)  have 
a  continuous  sample  derivative,  with  probability  one.  However,  this  assumption  was 
used  in  Part  A  of  the  proof,  but  not  at  all  in  Fart  B.  Hence  if  the  right  hand  side 
of  (1)  is  infinite,  the  equation  is  true  with  both  sides  infinite.  We  now  give  an 
example  of  a  case  where  the  integral  on  the  right  of  (1)  is  infinite,  and  hence  the 
corresponding  moment  is  infinite. 

For  this  example  we  take  a  covariance  function  of  the  form 

(16)  r (t)  =  1  -  X2  t2/2  -  x2/log | t |  +  o (t 2 / 1 og  | t | ) . 

That  this  can  be  done  follows  from  a  result  of  Pitman  [5].  In  fact  if 
H(X)  =  l-F(X)  +  F(-X)  for  X  >  0  we  can  choose  H(X)  so  that 

H(X)  ~  K/(X2  log2X)  as  X  — *>  - 
to  give  the  desired  form  (16). 

Consider  now  the  case  k  =  2,  and  u  =  0.  Then  one  can  show  by  some  calculation 

that 

**  — 

II  yly2  Pt(o,i)  dy1dy2  ~  K|Ar/(l-r2(x)) 

o  — 
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where  K  again  denotes  some  constant  and  A  is  the  covariance  matrix  for 

^(t^),  P  ( t  ^ P  * ( t  £ ) »  t  -  t?  -  t^.  But  straightforward  calculation  shows 

that 

|  A |  ~  t ^ / 1  og ^  | t  |  as  t  - •>  0 

and  hence 

m 

ff  yty2  Pt(°»i)  dy1dy2  ~  k/(M  log  | t I )  as  t  — >  o. 
o  - 

It  follows  from  this  that  the  right  hand  side  of  (1)  is  infinite,  in  this  case. 

Finally  we  note  a  sufficient  condition  for  P ( t )  to  possess  a  continuous  sample 
derivative,  with  probability  one,  is  tlat 

r(i)  =  1  -  X2t2/2  +  0(x2/ | log  It  !  |a  } 

for  some  a>  1.  This  follows  from  the  work  of  Belaev  [1].  In  our  case  r(t)  given 
by  (16)  just  fails  to  satisfy  this  requirement.  Hence  it  appears  that  the  require¬ 
ments  that  f  have  a  continuous  derivative  and  that  the  right  hand  side  of  (1)  be 
finite,  which  are  suf fi :ient  for  to  be  finite  and  given  by  (1),  are  also  very 
close  to  being  necessary  for  this  to  be  the  case. 
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Appendix. 

It  was  stated,  in  writing  down  certain  density  Junctions  that  if  t^  ..t^  are 
distinct  time  points,  then 

(i)  the  joint  distribution  of  £(t^) . . . f ( t^)  e ' ( t ^  )  .  .  .  £  '  ( t^)  is  non  singular,  and 

(ii)  the  joint  distribution  of  P(t^)  .  .  .1;  (t^)  is  non  singular. 

We  shall  now  prove  (i),  and  hence  (ii)  will  also  follow. 

Let  (as  assumed  throughout)  F(X)  have  an  absolutely  continuous  component  and 
write  A  =  f  ^  ]  for  the  covariance  matrix  of  ( t  ^ ) .  .  .  1  ( t^) ,  :  1  ( t  ^  )  .  .  .  c  1  ( t^)  .  Let 
A  ■  TAjj]  denote  the  covariance  matrix  of  ‘  ( t  ^  )  .  .  . !  ( t^)  B  =  [B^  ]  that  for 

s ' (t.  )  .  .  .  P ' (t  ) ,  and  C  the  matrix  of  "cross"  covariances,  C  .  =  cov(E(t.),  :’(t.)). 
Ik  i  J  1  j 

Then 


Let  O'  denote  the  vector  [^,...0,,  . 

—  1  k  ’  I 

which  are  not  all  zero.  Then  we  have 


$^1,  where  PF  ,  are  complex  numbers 


Au  ‘ 1  e 


i(tj-tj,)X 


dF(X) 


'j  l 


i(t  -t.)X 

■/  iX  e  J  dF(X) 


9  Kt  -t -)X 

Bjjf  =  /  X  e  J  '  dF(X) 


From  this  we  see  that 

it  X  2  it  X  2 

O' A 0  -  /  f  |  TJP  e  J  |  +  X  |.~*  e  J  | 

j  J  j  J 


-  iXSO.e  JZ$.e 

j  J  l  J 


*”i^t.  t 

+  iXSfcP.e  JED.e  ]  d  F  (> ) , 
;  j  t  < 


in  which  a  *  denotes  complex  conjugate.  Thus 

k  it , X  k  i t , X  2 


6'AB  m  f  |  1  0  e  j  +  iX  7  t .  e  j  |  dF(X) 

J-l  J  J-l  J 


11 


Now  since  the  t,  are  distinct  and  0.,  <t ,  not  all  zero  it  follows  that 
J  J  J 

k  it.X  k  it.X 

0 ,  e  ^  and  -i  X  "  $ .  e  ^ 


j  =  l  ’ 


i=l 


J 


are  different  regular  functions  of  the  complex 
variable  >,  and  hence  cannot  be  equal  for  more  than  a  countable  number  of  values  of 
X.  Hence  we  must  have  P'A  ~>  0  since  F(X)  has  an  absolutely  continuous  component 
and  the  measure  it  defines  is  not  concentrated  on  a  countable  set.  Thus  A  is  a 
strictly  positive  definite  matrix  and  the  distribution  thus  defined  is  non  singular. 

Finally  we  note  here  that  the  above  argument  can  be  easily  generalized  to 
include  an  arbitrary  number  of  derivatives.  That  is  if  F(\)  has  an  absolutely 
continuous  component  and  is  such  that  e(t)  has  n  sample  derivatives  ?(t)  | 1 (t ) . . . : ( t ) , 
then  for  any  distinct  tpt  ....t^,  the  joint  distribution  of 


=(n)(tl)...=(n) 


(t^)  is  non  singular. 
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